Boundary and finite-size effects in the competition between indirect magnetic 

exchange and Kondo screening 



(N 

o 

> 
O 

(N 



u 

-f— > 
-4-* 



I 

C 

O 

o 



> 
o 



(N 



Irakli TitvinidzeQ and Michael Potthoff 
I. Institut fur Theoretische Physik, Universitat Hamburg, Jungiusstrafie 9, 20355 Hamburg, Germany 

A system of conduction electrons can mediate an indirect magnetic exchange between magnetic 
impurites. The nonlocal exchange typically competes with the local Kondo screening of the impurity 
magnetic moments. In case of magnetic adatoms on non-magnetic surfaces this competition is 
expected to be affected by different confinement effects. Here we study this situation in the regime 
of intermediate hybridization strengths by means of real-space dynamical mean-field theory for a 
one-dimensional two-impurity Anderson model. Depending on the presence of a boundary and 
depending on the system size, a crossover between nonlocal and local singlet formation is observed 
in corresponding magnetic susceptibilities. The crossover is driven by the strong sensitivity of the 
Kondo temperature on the local density of states. 
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Introduction. Due to recent progress of experimen- 
tal techniques to investigate magnetic adatoms on sur- 
faces of non-magnetic metals not only the structural and 
electronic properties of adatoms and of the surrounding 
host material can be accessedjput also tailored magnetic 
model systems can be build p|. One of the most fasci- 
nating effects in this context is the interplay between an 
indirect magnetic exchange of the adatoms via the sub- 
strate electrons, i.e. the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction 0, on the one hand and Kondo 
screening Q of the adatom magnetic moment by the 
conduction-band electrons of the substrate on the other. 
Impurities, step edges and different confinement effects 
may result in strong variations in the local substrate den- 
sity of states. This strongly affects the Kondo screening 
as well as the RKKY interaction. Discrete eigenmodes 
can be found e.g. in quantum corrals located on a metal- 
lic surface |4]. A quantum box, i.e. a quantum system 
confined in all three dimensions, represents the extreme 
case. Here, the local density of states is singular and con- 
sists of a set of delta-functions. A quantum box can be 
realized, e.g., by semiconductor quantum dots [5|, or as 
shortened carbon nanotubes @, or as a single, e.g. metal- 
locene, molecule Q . Theoretically, the Kondo-box prob- 
lem has been studied extensively [§4lfl|. while there is 
not so much is known on the competition between Kondo 
screening and RKKY interaction in a quantum box [l2j . 

With the present contribution we report on studies of 
the effects of indirect magnetic exchange and its com- 
petition with Kondo screening for small one-dimensional 
model systems. We will consider systems with open and 
periodic boundary conditions as well as thermodynami- 
cally large systems to disentangle boundary and finite- 
size effects. Our main result is that there is a non- 
trivial competition between indirect magnetic exchange 
and Kondo screening which, at intermediate coupling 



strengths, is controlled by the local density of states 
rather than by finite-size gaps. 

Model. We consider the two-impurity Anderson model 
(TIAM): 
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Here /t a and cj a are creation operators of an elec- 
tron with spin projection a =t>-| a t the impurity (or 
"adatom") sites a — 1,2 or at the bath ("substrate") 
sites i = 1,2, ...,L, respectively. The corresponding 
occupation-number operators are n£ a = /t a f a a (n£ = 

"lt +n li) and n i° = c U c i,<T_( n i = n it+ n W- 1 is the 

hopping amplitude between neighboring substrate lattice 
sites and is used to fix the energy scale, i.e. t = 1. U and 
e denote the on-site Hubbard interaction and the local 
on-site energy for the adatom sites, [i is chemical poten- 
tial. V is the hybridization between an adatom site a 
and the nearest-neighboring substrate lattice site which 
is denoted by i a . We consider the particle- hole symmet- 
ric model with [i = and e = — U/2 at half- filling which 
implies ground-state expectation values (n£) = 1 and 
(nf) = 1 for all a and i. 

In the Kondo limit of the TIAM, where charge fluctua- 
tions on the adatom sites are strongly suppressed, there is 
an antiferromagnetic local exchange coupling J ~ V 2 /U 
between the adatom spin S Q and the substrate spin Si a 
by which one can characterize two relevant energy scales 
0, H| : Jrkky oc J 2 is the nonlocal indirect RKKY inter- 
action and Tk cx exp(— 1/| J|) the Kondo temperature. In 
the Kondo regime for Tk 3> | Jrkky| the local magnetic 
moments at the adatoms are independently screened by 
the substrate electron spins. On the other hand, in the 
RKKY limit for Tk <C | Jrkky | and after integrating out 
the substrate degrees of freedom, the problem reduces 
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to a two-spin Heisenberg model with exchange coupling 



^rkky = ^Xii™ , where Xi^° is the nonlocal static 
susceptibility of the substrate. 

In the present study we will restrict ourselves to the 
case of intermediate hybridization strength V to address 
a parameter regime where there are non-trivial effects of 
indirect magnetic exchange and which at the same time is 
clearly beyond the range of standard perturbative RKKY 
theory. To study systems of different size L, ranging 
from small quantum boxes of several tens of sites up to 
thermodynamically large systems where finite-size effects 
are absent, we employ the real-space dynamical mean- 
field theory (R-DMFT). 

Method. R-DMFT [l3[ is a comprehensive, thermo- 
dynamically consistent and non-perturbative mean-field 
approximation for correlated multi-impurity or lattice- 
fermion models with missing or reduced translational 
symmetry. The approach assumes the self-energy to be 
a local quantity but, in contrast to the standard DMFT 
[l4j, retains its site-dependence. Each inequivalent and 
correlated site is self-consistently mapped onto a single- 
impurity Anderson model. The set of Anderson models 
can be solved independently to get the local but site- 
dependent self-energies. The different sites are coupled 
via Dyson's equation for the multi-impurity model the 
solution of which provides us with the local Green's func- 
tions. These are used to determine the parameters of the 
impurity models in turn by means of the R-DMFT self- 
consistency equations. Details of the approach can be 
found in Refs. pJ,[H|. 

In a preceding study [TEj for small one-dimensional 
systems (L — 50) with open boundary conditions we 
have benchmarked the R-DMFT predictions against nu- 
merically exact density-matrix renormalization group 
(DMRG) calculations. The mean-field approach turned 
out to be even quantitatively reliable in the case of inter- 
mediate and strong hybridization V. 

For small V, however, it was found that the R-DMFT 
breaks down: Here the system becomes very susceptible 
to an externally applied magnetic field, and the mean- 
field approach artificially predicts a phase transition to a 
state with spontaneously broken symmetry. Specifically, 
R-DMFT yields a mixture of degenerate states with or- 
dered adatom local moments, | fj-i) an d | rather 
than the nonlocal singlet (| t,l) - I 4-,t))/V2 that is 
actually found in DMRG. 

Calculations. To study the TIAM at half-filling and 
zero temperature, the exact-diagonalization method [l6[ 
is used as an impurity solver. The impurity Green's func- 
tion is computed by means of the Lanczos technique for 
effective impurity models with n s = 10 sites. We have 
checked that bath-discretization errors are negligible for 
the results discussed below. To avoid a Kramers degener- 
acy of the ground state, we consider systems with an even 
number of substrate sites L. Calculations are done for 
one-dimensional systems with periodic boundary condi- 
tions. Small systems with L = 50 and L = 52 substrates 
sites as well as thermodynamically large systems are con- 
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FIG. 1: Local and nonlocal adatom susceptibilities xn and 
X21 as functions of V 2 /U for nearest-neighbor distance d = 1 
and U = 8. Results are shown for different system sizes as 
indicated. The results for L = 50 sites with open boundary 
condition are taken from Ref. [l5[ . Here the adatoms are 
located at positions symmetric to the chain center. Insets: 
results plotted on a logarithmic scale. 



sidercd. L = 10000 sites turned out to be sufficient to 
suppress finite-size effects in the calculation of the sub- 
strate Green's function (see Ref. [l5| for details). To dis- 
tinguish between finite-size and boundary effects we also 
compare the results with those of Ref. [I5[ where open 
boundary conditions have been employed. In that case 
the two adatoms are coupled to substrate sites located 
symmetric to chain center. 

Observables. To study the magnetic properties of 
the system, we consider the local adatom susceptibilities 
Xn = X22 and the nonlocal adatom-adatom susceptibili- 
ties X12 = X21 as weu as t ne adatom-substrate suscepti- 
bilities x-" b - These are obtained as numerical derivatives 
of the respective local magnetic moment with respect to 
a local magnetic field hp coupling the z-component of the 
local adatom spin to the system via % — > % — hpSp jZ : 



dml 



Xa/3 



dhp 



sub 

Ai/3 



dm c i 



h B =0 



dhp 



(2) 



h a =0 



Here m f a = {S a , z ) = |((nf )t ) - (ra^)) and ml = (s a , z ) 



)) are magnetic moments on the adatom 

1,2 L), 



(a = 1,2) and on the substrate sites (i 
respectively. 

Dependence on the coupling strength. We first consider 
the case where the two adatoms are coupled to neighbor- 
ing substrate sites, i.e. the distance between the adatoms 
is d = 1 (see Fig.[T]). In the regime with antiferromagnetic 
Kondo coupling V 2 /U > t, the susceptibilities for small 
systems (L — 50 and L = 52) are basically indistinguish- 
able from those for a thermodynamically large system 
(L = 10000). Contrary, for smaller values of V 2 /U, the 
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FIG. 2: Local and nonlocal adatom susceptibilities Xn an d X21 as functions of the distance d between adatoms for V 2 /U = 0.25 
and (7 = 8. Left: thermodynamically large system (L — 10000). The dashed line indicates the single-impurity susceptibility 
Xil = 6.54. Right: small systems (L = 50 and L = 52) for periodic and for open (taken from Ref. jig]) boundary conditions. 



results start to depend on the size of the system and on 
the type of boundary conditions. 

This trend can be explained by referring to the concept 
of the Kondo cloud. In the strong-coupling limit V 2 /U 
t, the size of the Kondo cloud £ is much smaller than size 
of the substrate £ <C L. Hence, finite-size effects in xn, 
which roughly is proportional to £, are negligible. For 
distance d — 1 boundary effects for open chains do not 
play a role either. 

A further decrease of V 2 /U finally leads to a divergence 
of the susceptibilities (see Fig. [1} . This must be seen as 
an artifact of the R-DMFT and should be interpreted as 
an indication for the crossover to the perturbative RKKY 
regime. As Fig. [1] shows, the position of this crossover is 
strongly affected by boundary and finite size effects. 

Here, we recall that besides the coupling strength, the 
Kondo temperature crucially depends on the local den- 
sity of states. For a finite system, where the local density 
of states is composed of delta-peaks only, the relevant 
quantity is the weight \Ui a k F \ 2 of the highest occupied 
one-particle energy eigenstate \kp) of the non-interacting 
substrate at i a , where \k) — L^ 1 / 2 J2iUik\i) and \i) is 
the basis orbital at site i. The weight is site-dependent 
in case of open chains. For L = 50 one can easily see 
that \Ui a k F 1 2 is larger compared to the weight for the in- 
finite system. Consequently, the Kondo temperature is 
higher than for the infinite system. This explains, why 
the crossover to the RKKY regime takes place for smaller 
values of V 2 /U in case of open boundaries and L = 50. 

For L — 52 sites and open boundaries (not shown in 
Fig.[TJ, we have just the opposite picture, and the Kondo 
temperature lower than for the infinite system. The dif- 
ference between L = An and L = An + 2 (for integer 
n) results from scattering at the chain boundaries caus- 
ing Friedel oscillations in the weight with a wave vector 
2kp = it. 



Contrary, for systems with periodic boundaries, the 
weight \Ui a k F \ 2 — 1 is site-independent and also inde- 
pendent of the system size. Still there is a clear oscilla- 
tion visible in the susceptibilities between systems with 
L = An and L = An + 2 sites. Here, the decisive differ- 
ence is that there is a (twofold degenerate) eigenenergy 
at zero excitation energy for L — An and thus a stronger 
effective coupling of the adatom moment and a higher 
Tk while there is a gap in the local density of states for 
L = An + 2. 

We conclude that, depending on the hybridization 
strength, finite-size and boundary effects show up in lo- 
cal adatom susceptibility xii which can be explained by a 
competition between the Kondo effect and indirect mag- 
netic exchange. As a divergency of xil also implies X21 
to be divergent, the same trends are also found there. 
Anyway, in the RKKY limit J — > 0, we have the sum 
rule xn + X21 = 0, and the behavior of X21 (see Fig. [T]) 
can also be interpreted as being a reminiscence to the 
RKKY limit. 

Infinite system: distance dependence. Fig. [2] (left) dis- 
plays our results for the local and the nonlocal adatom- 
adatom susceptibilities as a function of the distance d for 
fixed V 2 /U = 0.25. At this coupling strength, L = 10000 
substrate sites are sufficient to simulate the infinite sys- 
tem. Here, and in the following, we always consider odd 
distances d which results in an antiferromagnetic indirect 
exchange. The ground state is a spin singlet. 

The modulus of X21 decreases with increasing distance 
d, contrary to what would be predicted by an effec- 
tive RKKY two-spin model since a decreasing coupling 
(Jrkky fx 1/d in one dimension at half-filling) would 
lead to more susceptible moments. Here, rather a pic- 
ture of two separate Kondo clouds applies where with in- 
creasing d the adatom-substrate coupling becomes more 
important. Recall that for a gapped system we trivially 
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distance d 

FIG. 3: xn an< i X21 as functions of d for U = 8 and L — 
50 using periodic boundary conditions. Results for different 
V 2 /U as indicated. Note that the susceptibilities are rescaled 
with their value for d = 1. 



have the sum rule xn + X21 + J2i xfi b = 0. The large- 
d limit is equivalent with the single-impurity (SI) limit. 

This implies xn -> X11 and J2iXli h ~> Ei X™ b ' SI with 

Xu + Xa b ' SI = and therefore X21 -> 0. 

The local adatom susceptibility xn , on the other hand, 
increases with d and very slowly converges to Xn- We 
find X11 = 6-54 (see dashed line in Fig. [2] (left)) from a 
non-selfconsistent calculation for a single adatom. 

Finite systems: boundary effects. Fig. [5] (right) shows 
corresponding results for small systems. We first dis- 
cuss the results obtained for open boundaries (see blue 
lines). While the nonlocal adatom-adatom susceptibil- 
ity X21 again tends to zero with increasing d, we find 
a non-monotonic behavior of the local adatom suscep- 
tibility xn- As detailed in Ref. [lH, the decrease of 
X11 for large d is a boundary effect. With increasing 
d the adatoms move to the chain edges. This implies 
an increase of the local substrate density of states at the 
adatom sites which gives rise to an increase of the lo- 
cal Kondo temperature and thus a decrease of the local 
susceptibility. 

Finite systems with periodic boundaries. For periodic 
boundary conditions, the situation should be more simple 
as edge effects are absent. The results obtained for sys- 
tems with L = 52 substrate sites are indeed qualitatively 
similar to those for the infinite system (Fig. [5] (right), see 
brownish lines). Note that the maxima seen in xu and 
X21 at d — [L + l)/2 simply result from the equivalence 
between distances d and L — d. 

For L — 50, on the other hand, the nonlocal suscep- 
tibility X21 behaves quite different. Its modulus |x2i| 
increases with increasing distance for d < (L + l)/2 
(except for the extreme case of nearest neighbors, i.e. 
|X2i (d = 3)| < |X21 (d = 1)1)- An increase of |x2i| with 
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FIG. 4: Local adatom-substrate susceptibility Xnij where i\ 
is the site "below" the first adatom, as a function of the dis- 
tance d to the second adatom. Results for V 2 /U = 0.25 and 
U — 8. Upper panel: large system, L — 10000. Lower panel: 
small systems, L — 50 and L — 52, with periodic and open 
(taken from Ref. Q^) boundary conditions. 



increasing distance between the adatoms could be under- 
stood perturbatively within RKKY theory as mentioned 
above. However, V 2 /U = 0.25 and U — 8 would re- 
fer to the strong-coupling limit (J = 2 in units of t) 
and, anyway, there are still significant charge fluctua- 
tions at the adatom sites, and hence the RKKY picture 
does not apply. Nevertheless, it is tempting to interpret 
the trend of |x2i| as a reminiscence to the perturbative 
RKKY regime. Namely, as discussed above, the system 
has a lower Kondo temperature compared to the infinite 
system and compared to the case L = 52, too. That the 
formation of a nonlocal singlet is favored here, is also cor- 
roborated by the larger absolute values of X21 and by the 
fact that X21 behaves almost complementary to xn as it 
is characteristic for the RKKY limit where xn +X21 = 
holds. 

Crossover between local and nonlocal singlets in finite 
systems. As a further check we have performed cal- 
culations for L — 50 using periodic boundary condi- 
tions again but at larger values of V 2 /U (see Fig. [3] 
and note that the results are rescaled for better compari- 
son). In fact, with increasing V 2 /U the nonlocal adatom- 
adatom susceptibility X21 nicely shows a crossover from 
the regime where nonlocal singlet formation dominates 
to a Kondo regime with local singlets. This is indicated 
by the qualitative and systematic change of the d depen- 
dence of X21 when increasing V 2 /U. 

With this we can also understand the anomalous be- 
havior of the susceptibility at d = 1 compared to d = 3: 
Due to the proximity to the RKKY regime, the feed- 
back of nonlocal magnetic correlations on the selfenergy 
is increased. Since this is neglected within R-DMFT, the 
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behavior is most probably an artifact of the mean-field 
approximation. Consistent with this are the slight de- 
viations of the R-DMFT results from numerically exact 
DMRG data that have been observed in Ref. [l5[ at cl = 1 
for similar model parameters. 

Finally, Fig. |4] displays the results for the local adatom- 
substrate susceptibility xf"i which measures the response 
of the substrate site i\ "below" the adatom to a local field 
at the same adatom. This probes the strength of local 
Kondo singlet formation, xt^i m f ac t is negative in all 
cases. Consistent with the above interpretation the mod- 
ulus |xf"i1 is the largest for the system with the lowest 
7k, i.e. for L — 50 with periodic boundaries. For open 
boundaries, we again have a site-dependent Kondo tem- 
perature. This explains the non-monotonic behavior of 
|x|"i | as a function of the distance to the second adatom 
(note that for periodic boundaries we trivially have an 
equivalence between distances d and L — d). We also 
find the artifacts of the mean-field approximation to be 
the strongest for systems with dominant nonlocal singlet 
formation, i.e. for L = 50 with periodic boundaries. 

Conclusions. We have studied model systems where 
two correlated sites with well-defined local magnetic mo- 
ments interact locally with the itinerant moments of non- 
interacting conduction electrons. As has been verified 
previously by comparison with DMRG calculations, a 
regime of intermediate hybridization strengths is well ac- 
cessible to real-space dynamical mean-field theory. Here 
we find that just in this parameter range a non-trivial 
competition between nonlocal indirect exchange and lo- 



cal Kondo screening takes place which is strongly affected 
by finite-size and boundary effects. This is most clearly 
seen in the static nonlocal inter-impurity magnetic sus- 
ceptibility. Depending on the type of boundary condi- 
tions and depending on the system size, a crossover be- 
tween different distance dependencies that are character- 
istic for the two extreme cases of the nonlocal RKKY 
singlet and local Kondo singlets respectively, is found as 
a function of the coupling strength. Rather than a low- 
energy cutoff due to finite-size gaps, the main mechanism 
controlling this competition at intermediate hybridiza- 
tion strength is the exponentially strong dependence of 
the Kondo temperature on the local density of states. 
The LDOS in turn is sensitively affected by the system 
size and according odd-even effects and, via Friedel oscil- 
lations, by the presence of the system's boundary. While 
a one-dimensional model system has been the focus of 
the present study, real-space DMFT can easily be applied 
to basically arbitrary structures of magnetic adatoms on 
surfaces of higher-dimensional lattices. Corresponding 
studies are intended for the future. 
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